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Abstract. This paper proposes a new class of copulas which character- 
ize the set of all twice continuously differentiable copulas. We show that 
our proposed new class of copulas is a new generalized copula family that 
include not only asymmetric copulas but also all smooth copula families 
available in the current literature. Spearman's rho and Kendall's tau 
for our new Fourier copulas which are asymmetric are introduced. Fur- 
thermore, an approximation method is discussed in order to optimize 
Spearman's rho and the corresponding Kendall's tau. 



1. Introduction 

Recently, a study of dependence by using copulas has been getting more 
attention in the areas of finance, actuarial science, biomedical studies and 
engineering because a copula function does not require a normal distribu- 
tion and independent, identical distribution assumptions. Furthermore, the 
invariance property of copula has been attractive in the finance area. But 
most copulas including Archimedean copula family are symmetric functions 
so that the copula model fitting for asymmetric data is not appropriate. 
Liebscher (jEj) introduced two methods for the construction of asymmetric 
multivariate copulas. The first is connected with products of copulas. The 
second approach generalizes the Archimedean copulas. The resulting cop- 
ulas are asymmetric but are little extension of the parametric families of 
copulas. This paper proposes a new generalized copula family which include 
asymmetric copulas in addition to all copula families available in the current 
literature. 

We will characterize the set of all twice continuously differentiable copulas 
that may arise. We will start with some basic concepts of copulas in the 
next section. In Section 3 we characterize a class of differentiable copulas 
and state our main theorem. In Section 4 we will discuss some well-known 
copulas and introduce a new class of copulas (asymmetric in general) in sup- 
port of our main result. Finally, in Section 5, we will study the dependence 
structure by calculating Spearman's rho and Kendall's tau and describe a 
method to maximize and minimize Spearman's rho using an approximation 
technique for a special class of copulas that arises from our construction. 
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2. Definitions and Preliminaries 

In this section we recall some concepts and results that are necessary to 
understand a (bivariate) copula. Throughout this paper I denotes the unit 
interval [0, 1]. A copula is a multivariate distribution function defined on I n , 
with uniformly distributed marginals. In this paper, we focus on bivariate 
(two-dimensional, n = 2) copulas. 

Definition 2.1. A bivariate copula is a function C : I 2 — > I, which satisfies 
the following properties 

(PI) C(0,v) = C(u,0) =0, Vu,v£l 

(P2) C(l,u) = C(u,l) = u, VuGl 

(P3) C(u 2 ,V 2 ) + C(ui,Vl) — C(ui,V 2 ) — C(U2,V\) > 0, VUi, U 2 ,V\,V 2 G 

I with u\ < u 2 , v\ <V2- 

The importance of copulas has been growing because of their applica- 
tions in several fields of research. Their relevance primarily comes from the 
following theorem of Sklar (see (Sj): 

Theorem 2.2. Let X and Y be two random variables with joint distribution 
function H and marginal distribution functions F and G, respectively. Then 
there exists a copula C such that 

H(x,y) = C(F(x), G(y)) 
for all x,y G M. If F and G are continuous, then C is unique. Otherwise, 
the copula C is uniquely determined on Ran{F) x Ran(G). Conversely, if 
C is a copula and F and G are distribution functions, then the function H 
defined above is a joint distribution function with margins F and G. 

Sklar's theorem clarifies the role that copulas play in the relationship 
between multivariate distribution functions and their univariate margins. A 
proof of this theorem can be found in |SS| . 

Definition 2.3. Frechet lower and upper bounds for copulas are denoted by 
Cl and Cjj, respectively, and defined by 

Cl(u, v) = max{-u + v — 1, 0}, 

Cu(u,v) = min{n,ti}, 

for all (u, v) G I 2 . 

Then it well-known that for any copula C(u,v), 

C L {u, v) < C(u, v) < Cu(u, v), V(«, v) G I 2 . 

In this paper we will mainly concentrate on copulas which are twice con- 
tinuously differentiable, i.e., C(-,-) G C 2 (I 2 ). With this assumption and 
using property 3 (P3) of copulas, for all u\, u 2 , v\, V2 G I with u\ < u 2 , v± < 
i>2, we have C{u2,V2) — C(ui,v 2 ) > C(u 2 ,vi) — C(ui,v\). This implies 
C(u 2 , •) — C(u\, •) is monotonically increasing in the second variable. Hence 
^ [C(u 2 , v ) — C{u\, v)] > 0. Therefore ^C(-,v) is increasing in its first 
variable. Hence we deduce the following lemma: 
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Lemma 2.4. The following two statements are equivalent: 

i. C is a twice continuously differentiable copula. 

ii. C satisfies the following Dirichlet inequality problem: 

with boundary conditions: 

C(u,0) = C(0,v) = 0, 

C(u,l) = C(l,u) = u. 

3. Characterization of C 2 copulas 

In this section we will start by solving the above problem stated in Lemma 
|2,4| and then we will characterize all twice differentiable copulas. 

Suppose 7 : I 2 — > M is a continuous real-valued function. Then inequality 
can be reformulated as follows, 

d 2 

q u q v C ( u ^ v ) = 7 2 ( u > v )- 

Integrating twice we get 



C{u,v) 



/ / 7 2 (s,t) ds dt + H(u) + G(v) 
Jo Jo 



where H and G are two arbitrarily real- valued functions of u and v, re- 
spectively. Using boundary conditions C(u,0) = C(0,v) = 0, we have 
H(u) = —G(v) = constant. Hence C has the following form 

(2) C(u,v) = l 2 (s,t) ds dt. 

Jo Jo 

Now using the boundary condition C(l, v) = v, we have 



7 (s, t) ds dt = v. 

'o 

Differentiating both sides with respect to v, we have 

(3) f -y 2 (u,v) du = 1, \/v G I. 

Jo 

Similarly using the fourth boundary condition C(u, 1) = u, we have 



(4) 

This leads to the following theorem 



/ j 2 (u,v) dv = 1, Vtt € I. 
Jo 
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Theorem 3.1. Suppose h : I 2 — > [— l,oo) is a continuous real-valued func- 
tion such that 

(5) f h(u, v) dv = V«G I, 
Jo 

(6) / h(u,v) du = to € I, 
J o 

(7) C(«,t;) = / / l + h(s,t) dsdt. 

Jo Jo 

is a copula. Furthermore, every twice continuously differentiable copula is 
of the form given in 

Proof. If we substitute 1 + h(u, v) for 7 2 (u, v) in Equations pj and Q, 
it is easy to verify that every twice continuously differentiable copula is given 
by 0. 

n q2 C 

To prove the other direction, we have from (7), - - = 1 + h(u,v), 

|_[ ouav 

which is continuous and non-negative on I 2 . It is easy to check C(u,Q) = 
C(0,v) = for all (u,v) G I 2 . We also have, 



pv pi 

C(l,v) = l + h(s,t)dsdt 

Jo Jo 



v 

1 dt 



o 



V. 



Similarly, we can show that C(u,l) = u. Hence by Lemma 2.4 C is a 



copula. □ 



Remark 3.2. One importance of Theorem 3.1 is the fact that in general 
there is no assumption of symmetry on h(u,v) and hence on C(u,v). With 
the help of above theorem, in the next section we will construct a class of 
examples of non-symmetric copulas. 

4. Examples 



It is quite easy to verify Theorem 3.1 for well-known C copulas by con- 



structing the corresponding h. In this section we will begin with showing 
those renowned examples and later we will show how to construct other 
copulas. 

4.1. Archimedean Copulas. Archimedean copula is a very interesting 
class of copulas, whose investigation arose in the context of associative func- 
tions and probabilistic metric spaces (see |SSj ) and today has also many 
applications in the statistical context (see [N]). Moreover, Archimedean 
copulas are widely used in applications, especially in finance, insurance and 
actuarial science, due to their simple form and nice properties. 
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Definition 4.1. Let </?:!—>• [0, oo] be a continuous, decreasing function 
such that = 0. The pseudo-inverse of cp is the function denoted i/jH 
with domain [0, oo] , range I and defined by 

^ 1 J ~ \ i/ y>(0) < t < oo 

Definition 4.2. Lei <£> : 1 — > [0, oo] be a continuous, convex, strictly decreas- 
ing function such that = 0. Then a copula C is called Archimedean if 
it can be written as 

C{u, v) = ip^ (<p(u) + tp(v)) , V(u, v) G I 2 . 

ip is called an additive generator of C . 

Now assuming C G C 2 (I 2 ), let us define 

h(u, v) = (/j' -1 ! (ip(u) + tp(v)) (p'(u)<p'(v) — 1. 
Then it is evident that h is continuous. Also one can easily show that 

f'U rv 

C(u,v) = / 1 + h(s,t) dt ds and h satisfies both (5) and (6). 

Jo Jo Li L 

Example 4.3 (Frank Copulas). The Frank copula is an Archimedean copula 
given by 

C(u, v) = -- In |l + * ' 

where 9 G K \ {0} is a parameter. Its generator is given by 

( £-® x — i " 

V? e (x) = -ln 

T/ien 



1 



;i_ e e + e 0+^)2' ™«v e - te -l' 
Therefore we have 

_ 9e^ l+u+v \e e - 1) 



One can easily verify that h satisfies Theorem 3.1 



4.2. FGM Copulas. The Farlie-Gumbel-Morgenstern (FGM) copula takes 
the form 

C(u,v) =uv(l + 6(l-u)(l-v)), 

where 9 G [—1, 1] is a parameter. The FGM copula was first proposed by 
Morgenstern (1956). The FGM copula is a perturbation of the product 
copula; if the dependence parameter 9 equals zero, then the FGM copula 
collapses to independence. This is attractive primarily because of its simple 
analytical form. FGM distributions have been widely used in modeling, for 
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tests of association, and in studying the efficiency of nonparametric pro- 
cedures. However, it is restrictive because this copula is only useful when 
dependence between the two marginals is modest in magnitude. 

Let h(u,v) = 6(1 — 2u)(l — 2v). It can be easily verified that h agrees 
with Theorem 3.1 and it generates the FGM copulas. 

4.3. Fourier Copulas. The following lemma will introduce a new class of 
C 2 copulas. 

Lemma 4.4. Suppose (a n ) , (b n ) , (c n ) , (d n ) G (r , the space of sequences 
whose series is absolutely convergent, are sequences of real numbers. Also 
suppose h is a real-valued function on I 2 defined by 

oo oo 

h(u, v) = ( a « cos (2irnu) + b n sin (2irnu)) (c m cos (27rmv) + d m sin (2Timv)) , 
n=l m=l 

then 

l 

h(u, v) du = 0, \/v G I, 

l 

h(u, v) dv = 0, Vn G I. 

o 

oo 

Furthermore, if \/ a n + b„ \J c 2 n + d^ < 1, then h(u,v) > —1 for all 

n,m=l 

(u,v) G I 2 . 

Proof. Notice that we can rewrite h(u, v) as, 

h(u, v) = a nCm cos (2imu) cos (2-Kmv) + a nd m cos (2-irnu) sin (2irmv) 

n,m n,m 

+ b n c m sin (27rnu) cos (2irmv) + b n d m sin (2imu) sin (2irmv). 

n,m n,m 

The first conclusion follows from the fact that the sequences (a n ), (b n ), (c n ), (d n ) G 

i 1 and / sin (2imx) dx = cos (2irnx) dx = 0, Vn G N. 
Jo Jo 
Now to prove h(u, v) > — 1 for all (u, v) G I 2 , first notice that \/a 2 + 6 2 < 

n 

^ (|On| + |&n|) < OO. Also, 



~V a n J i~b 2 l < a n cos (2-kuu) + 6 n sin (27rnu) < -y/a^ + b 2 , 
for all u G I, n G N. Hence we have 

- Va-l + bWcm + dlt < h(u, v) < ^ \/a 2 + 6 2 ^c 2 , + d 2 ,. 

n,m n,m 

Therefore the additional hypothesis on h guarantees that h(u,v) > — 1. □ 
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It is evident that h(-,-) is continuous and therefore by Theorem 3.1, C, 
defined by 



C(u,v)= / / l + h(s,t) dsdt, 
Jo Jo 



is a copula, where h is of the form given in Lemma 4.4 
Noting that 

oo 

a n cos(27rnn) + b n sin(27rnn) = ^ 7ne 2?m "\ 

n=l n.eZ\{0} 

where a n = 7„ + 7_ n and b n = i{^ n - j- n ), or equivalently, 

Q"n ibn 

In = 7-n = ^ ) 

it follows that |7„| = ^ a \ +b " . Similarly, if S m = 6~ - ^4 



m — 2 



h(u,v)= InSm exp(27rz(nu + mti)), 

n,m£Z\{0} 

with sequences (7 n ), (<5 m ) G and ||7n||||^m|| < j> then h will satisfy 
the conclusions of Lemma 4.4, and will generate a copula by eq. ([7]). This 
restatement of Lemma 4.4 clearly indicates that h is obtained from products 
of functions in the disc algebra with vanishing zero-th moment and where 
the product of the i l norms of the coefficients are < \. More precisely, 

Theorem 4.5. If h is a function on the 2-torus arising from the product 
of functions on the unit circle, each section of h has a vanishing zero-th 
moment, and the product of i 1 norms of Fourier coefficients of components 
of h is < j, then h generates a C 2 -copula. 

This theorem provides a large class of examples from which one may 
construct copulas with optimal properties (p and r). 

Example 4.6. Let h\ = c\ = 1, b n = c n = 0, Vn ^ 1; a n = d n = Vn ; then 

h(u, v) = sin (2iru) cos (2irv). 

Define C as follows 

C(u,v) = / / 1 + sin (27rs) cos (2irt) ds dt 
Jo Jo 
[ v 1 

= / u-\ {1 — cos (27ru)} cos (2nt) ds dt 

Jo 27r 

= uv H ^{1 — cos (27rn)} sin (2irv) 

Then it is easy to verify that C forms a copula and more importantly it is 
not symmetric. A contour plot of C is given in Figure [7} 
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Figure 1. Contour plot of a Fourier copula in Ex. 4.6 



5. Spearman rho and Kendall tau 



The two most commonly used nonparametric measures of association for 
two random variables are Spearman rho (p) and Kendall tau (r). In general 
they measure different aspects of the dependence structures and hence for 
many joint distributions these two measures have different values. 

Definition 5.1. Suppose X and Y are two random variables with marginal 
distribution functions F and G, respectively. Then Spearman p is the ordi- 
nary (Pearson) correlation coefficient of the transformed random variables 
F(X) and G(Y), while Kendall r is the difference between the probability of 
concordance Pr[(Xl — X2)(Y1 — Y2) > 0] and the probability of discordance 
Pr[(Xl - X2)(Y1 - Y2) < 0] for two independent pairs (XI, Yl) and (X2, 
Y2) of observations drawn from the distribution. 

In terms of dependence properties, Spearman p is a measure of average 
quadrant dependence, while Kendall r is a measure of average likelihood 
ratio dependence (see [Nj for details). If X and Y are two continuous random 
variables with copula C, then Kendall r and Spearman p of X and Y are 
given by, 





(9) 
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5.1. p, t for Fourier copulas. From Lemma ( |4.4[ ), if we define 

h(u, v) = ( a n cos (2nnu) + b n sin (2imu)) (c m cos {2irmv) + d m sin (27rmv)) , 

n m 

where \/ a n + V c m + ^ 1) then the Fourier copulas are given by 



C(u,v) = 1 + h(s, t) ds dt 

Jo Jo 

= UV+} - — k [a n sin (2irnu) + b n {l — cos (27rnu)}] 

4:TT z nm 

n,m 

[c m sin {2irmu) + d m {l — cos (2irmv)}] 
Using ^ and Q, we have 

3 b, ,<],,, 



ir 2 '■ — ' nm 



2 \ b n d m 
7r 2 ' nm 



Now since - V ^ + ^^ + 4 < V] - 5 — < V] ^ + ^^ + 4. 

n.m 

we can conclude that 



nm 

n,m n,m 



-0.304 — =• < p < 4r ss 0.304, 



2 — r — 2 



-0.203 » = < r < rj 0.203. 

7T^ 7T^ 



Remark 5.2. Fourier copulas can be generalized by writing /i as follows, 
h(s,t) = ^2 a n,m exp(2iri(ns + mt)), 

n,m£Z\{0} 

where c^n,m = a —n,—m Vti, ?Tt £ Z \ {0} and ^ ^ |c*n,m| + |°!— n,m 

|< |. The 

n,m£N 

latter condition here is to ensure that h will have range in [—1, oo). Notice 
that for all n,m / 0, ot.- n ,m is equal to a n - m and hence no additional 
conditions are necessary to assure h to be real-valued. 

™ . . , , 3 \ -> ttn,m 2 \ -> OL n m 

1ms yields p = ^ > , and r = = > , and it 

n,m£Z\{0} n,m£Z\{0} 

can be shown again that —\ < p < \, and — \ < r < \. 
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5.2. Optimizing p for C 2 -copulas. In this section we will optimize Spear- 
man's rho using an approximation method for C 2 -copulas with h of the form, 
h(x,y) = (p(x)ifi(y), where ip and ip both are continuous real-valued func- 
tions on I. Notice that in this special case, p can be simplified into the 
following form, 



12 



P 



1 + ip(s)ijj(t)ds dt 

lt=0Js=0 
rl pu rl rv 

12 / tp(s)dsdu l I tp(t)dtdv 
Jo Jo Jo Jo 



du dv 



This suggests that optimizing p is equivalent to optimizing both L J* " (p(s)ds du 
and Jq Jq ip(t)dt dv. 

Define G(u) := (p(s)ds and H{v) := f£ ip(t)dt. Then for some positive 
o-i, «2, A, 02, the optimization problems become, 
■ i 



max/min 
subject to 



h 



G(u) du 



G(0) = G{1) = 
-«i <G'(u) <h, 



max/min 
subject to 



H(v) du 



H(0) = H(l) 

— QJ2 < H (v) 



= 
</3 2 - 



Although it apparently looks like these optimization problems can be 
solved independently, they are related by the fact that G'(u)H'(v) = ip(u)vjj(v) 
— 1. This implies min{— ai/32, — «2/3i} > —1- For the optimal possibility, we 
choose, /?2 = — (ai) _1 and ai = — (/3i) . This is evident from the fact that 
both 1\ and I2 can be positive or negative, p max will occur either if both I\ 
and I2 are maximum or if both are minimum and p m i n will occur if one of 
I\ and I2 is maximum and the other is minimum. 

Geometrically, I± will be maximum if G has the form as in Figure [2] and 
will be minimum if G has the form as in Figure [3} 



> 




Slope 



Slope -a\ 



(1.2 0.4 

Figure 2. G, Maxi- 
mizing I\ 



Slope -a 



Slope /3, 




Figure 3. G, Mini- 
mizing I2 



One can easily prove that, in order to optimize I±, /3\ must be equal to a\. 
For convenience, now onwards we will write a for ot\. This suggests that if 
G{x) = GM(x) = -a\x - 0.5| + 0.5a, or G(x) = Gm(x) = a\x - 0.5| - 0.5a 
then I\ will be maximum or minimum, respectively. But in either case, G is 
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not differentiable at x = 0.5, and hence (p is not continuous. To avoid this, 
we will approximate G by a smooth function as follows: for arbitrarily small 
e > 0, define 

GM(x) = -Gm{x) = £ (Vl + 4e 2 - ^(1 - 2x) 2 + 4e 2 ) . 



as 



It is worth noting that sup \\GM(x) - GM(x)\,\Gm(x) - Gm(x)\\ -> 
e -)■ and -q < GM (x),Gm (x) < a. 




0.2 0.4 0.6 0.8 




0.2 0.4 0.6 0.8 



Figure 4. GM, GM 
for a = 5, e = 0.1 



Figure 5. GM, GM 
for a = 5, e = 0.03 



We can similarly optimize I2 by approximating maximum and minimum 
of H by the following functions 

HM(x) = -Hm{x) = — (Vl + 4e 2 - y^ 1 - 2 ^) 2 + 4e 2 ' 

Hence optimum values of /) will be obtained by approximating h by the 
following functions, 

•/ — / •/ •/ (1 - 2x)(l - 2v) 

/4ax(z, V) = GM (x)HM (y) = GM (x)HM (y) - " -' 



^(1 - 2x) 2 + 4e 2 V(l - 2y) 2 + 4e 2 ' 

(l-2x)(l-2y) 
" V(l - 2x) 2 + 4^(1 - 2y) 2 + 4e 2 ' 
Notice that each of /i^ ax and h £ { will generate a copula as it satisfies 



h £ rmn(x,y) = GM {x)Hm (y) = Gm (x)HM (y) 



all the hypothesis of Theorem 3.1 Then the corresponding Spearman's rho 
and Kendall's tau are given by, 



Pmax 



Pram 



l + 4e 2 -4e 2 coth _1 (\/l + 4e 2 ) 
l + 4e 2 - 4e 2 coth _1 ( + 4e 2 )) 2 



1 + 4e 2 + 4e 2 



1 + 4e 2 + 4e 2 



1 + 4e 2 - 2e 2 coth _1 (\/l + 4e 2 )J In ^ 
1 + 4e 2 - 2e 2 coth -1 (\/l + 4e 2 )) In 



l + 2e 2 - Vl + 4e 2 \ 
2^ J 

/ 1 + 2e 2 - Vl + 4e 2 
2^ 
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e 


Pmax 




T 6 
max 


mm 


1 


0.0726437 


-0.0726437 


0.0484292 


-0.0484292 


0.1 


0.644923 


-0.644923 


0.429949 


-0.429949 


0.01 


0.747539 


-0.747539 


0.498359 


-0.498359 


0.001 


0.749962 


-0.749962 


0.499974 


-0.499974 


0.0001 


0.749999 


-0.749999 


0.5 


-0.5 



Table 1. p and r values as e changes. 



The optimal values of p and corresponding r will be obtained by letting 
e — > 0. Table [T] shows how the values of p approach the optimal values as 
e and it is clear that -0.75 < p < 0.75 and -0.5 < r < 0.5. 

6. Conclusion 

We proposed a new generalized copula family which include not only 
asymmetric copulas but also all copula families available in the current lit- 
erature. Especially, the family of Fourier copulas we proposed is very useful 
copula family for analyzing asymmetric data such as financial return data 
or cancer data in Bioinformatics. In our future study, we will extend our 
copula method to a multivariate case and then incorporate time varying 
component to our proposed method. 
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